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Abstract 

We present some Pfaffian identities, which are completely different from the 
Pliicker relations. As consequences we obtain a quadratic identity for the number 
of perfect matchings of plane graphs, which has a simpler form than the formula by 
Yan et al (Graphical condensation of plane graphs: a combinatorial approach, The- 
oret. Comput. Sci., to appear), and we also obtain some new determinant identities. 
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1 Introduction 

Let A = ( 

0'ij)nxn be a skew symmetric matrix of order n and n is even. Suppose that 
71 = {{si, ti), {s2, ^2), • • • , , ^a)} is a partition of [n], that is, [n] = {si, ti} U {s2, ^2} U 
. . . U {sh, tn}, where [n] = {1, 2, . . . , n}. Define: 



2 

'''Sltl ) 



K = sgn{sitiS2t2 . . . J]^ 

1=1 

where sgn{sitiS2t2 ■ ■ ■ s^t^) denotes the sign of the permutation sitiS2t2 • • • s^tii. Note 
that depends neither on the order in which the classes of the partition are hsted nor 
on the order of the two elements of a class. So indeed depends only on the choice of 
the partition vr. The Pfaffian of A, denoted by Pf{A), is defined as 

Pf{A) = J2b., 
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where the summation is over all partitions of [n], which are of the form of vr. For the sake 
of convenience, we define the Pfaffian of A to be zero if A is a skew symmetric matrix of 
odd order. The following result is well known: 

Proposition 1.1 (Cayley Theorem, [Ij) For any skew symmetric matrix A = {aij)nxn 
of order n, we have 

det(A) = [Pf{A)f. 

Suppose that G = iy{G), E{G)) is a weighted graph with the vertex set V{G) = 
{l,2,...,n}, the edge set E{G) = {ci, 62, . . . , and the edge— weight function u : 
E{G) — > 71, where Lj{e) := Ue = dij (7^ 0) if e = (z, j) is an edge of G and Ue = = 
otherwise, and TZ is the set of real numbers. Suppose is an orientation of G. Let 
A{G^) = {bij)nxn be the matrix of order n defined as follows: 

{ttij if is an arc in G^, 
—aij if (j, i) is an arc in G^, 
otherwise. 

A{G^) is called the skew adjacency matrix of G^ (see fZ]). Obviously, A{G^) is a skew 
symmetric matrix, that is, {A{G^))'^ = — A(G^). 

Given a skew symmetric matrix A = {aij)nxn with n even, let G = {V{G), E{G)) be 
a weighted graph with the vertex set V{G) = {1,2,..., n}, where e = (z, j) is an edge of 
G if and only if aij 7^ 0, and the edge— weight function is defined as Ue = \aij\ if e = {i,j) 
is an edge of G and = otherwise. Define as the orientation of G in which the 
direction of every edge e = of G is from vertices i to j if a^j > and from vertices 
j to i otherwise. We call G^ to be the corresponding directed graph of A. Obviously, 
A = y4(G^). It is not difficult to see that the Pfaffian Pf{A) of A can be defined as 

Pf{A)= J2 ^- 

weMiG) 

where the summation is over all perfect matchings vr = {(si,)f:i), (52,^2), • • • , {s^,t!i)} of 
G, and b-,, is the product of all ^^(s,,*;) for 1 < i < |. 

Pfaffians have been studied for almost two hundred years (see ^1I2H] for a history), 
and continue to find numerous applications, for example in matching theory ^| and in the 
enumeration of plane partitions [2H]- It is interesting to extend Leclerc's combinatorics 
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of relations for determinants to the analogous rules for Pfaffians. By tools from 
multilinear algebra Dress and Wenzel |^ gave an elegant proof of an identity concerning 
pfaffians of skew symmetric matrices, which yields the Grassmann-Pliicker identities (for 
more details see Sect. 7). Okada pS] presented a Pfaffian identity involving elliptic 
functions, whose rational limit gives a generalization of Schur's Pfaffian identity. Knuth 
[T^ used a combinatorial method to give an elegant proof of a classical Pfaffian identity 
found in [21] ■ Hamel [H] followed Knuth's approach and introduced other combinatorial 
methods to prove a host of Pfaffian identities from physics in [7| |^ EU]. Hamel also 
provided a combinatorial proof of a result in [2ZI and a new vector-based Pfaffian identity 
and gave an application to the theory of symmetric functions by proving an identity for 
Schur Q-functions. For some related recent results see also jHl El UH 12^] • 

This paper is inspired by two results, one of which is that we can use the Pfaffian 
method to enumerate perfect matchings of plane graphs (see fU^]). Inspired by the 
Dodgson's Determinant-Evaluation Rule in |i] and the Pliicker relations for Pfaffians, 
Propp [21] , Kuo Pl] and Yan et al [HS] obtained a method of graphical vertex-condensation 
for enumerating perfect matchings of plane bipartite graphs. The second is that by using 
the Matching Factorization Theorem in Yan et al |S2] found a method of graphical edge- 
condensation for counting perfect matchings of plane graphs. It is natural to ask whether 
there exist some Pfaffian identities completely different from the Pliicker relations, which 
can result in some formulas for the method of graphical edge-condensation for enumerating 
perfect matchings of plane graphs. The results in Section 3 answer this question in the 
affirmative. As applications, we obtain two new determinant identities in Section 4.1 and 
we prove a quadratic relation for the number of perfect matchings of plane graphs in 
Section 4.2, which has a simpler form than the formula in |32j . 

2 Some Lemmas 

In order to present the following lemmas, we need to introduce some notation and 
terminology. If / is a subset of [n] , we use Aj to denote the minor of A by deleting rows 
and columns indexed by I. If / = {ii,i2, ■ ■ ■ ^ [n] and ii < i2 < ■ ■ ■ < ii, we use 
-P/a(^i^2 ■ ■ ■'ii) ='■ PfA{I) to denote the Pfaffian of ^[n]\/- Following Knuth's notation in 
|13j . for two words a and (3 we define s(a,/5) to be zero if either a oi (3 has a repeated 
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letter, or if (3 contains a letter not in a. Or, if these are not the case, s{a,(3) denotes 
the sign of the permutation that takes a into the word (where a\(3 denotes the 

word that remains when the elements of (3 are removed from a). Let S" be a subset of 
{1, 2, . . . , n}. We call S an even subset if IS"! is even and an odd one otherwise. 

Dress et al jS] used tools from multilinear algebra to prove a Pfaffian identity, which 
was found by Wenzel jHI], as follows: 

Lemma 2.1 (Wenzel [31] and Dress et al [3]) For any two subsets h, h ^ [n] of 
odd cardinality and elements ii,i2, . . . ,it € [n] with ii < 12. . . < it and {zi, i2, ■ ■ ■ , it} = 
/1A/2 =: {hXh) U (/2\/i), if A = {aij)nxn is a skew symmetric matrix with n even, then 

E(-1)^P/a(/iA{v})P/(J2A{v}) = 0. 

T=l 

A direct result of Lemma 2.1 is the following lemma, which will play an important 
role in the proofs of our main results. 

Lemma 2.2 Suppose that A = {a.ij)nxn is a skew symmetric matrix with n even and a is 
an even subset of [n]. Let jd = {ii,i2, ■ ■ ■ ,"^2?} ^ where ii < 12 < ■ ■ ■ < i2p- Then, 

for any fixed s G [2p], we have 

2p 

PfA{a)PfA{a(3) = Z{-iy^'^'PfA{msti)PfA{a(3\tsti), 
1=1 

where PfAictish) = 0. 

The following result is a special case of Lemma 2.2. 

Corollary 2.1 Suppose that A = {aij)nxn is a skew symmetric matrix and {i,j,k,l} C 
[n] . Then 

P/(A|„-,,,})P/(A) = P/(A{,,|)P/(%,,|)-P/(A|,,}) 

(2.1) 

Remark 2.1 There exists a similar formula on the determinant to Corollary 2.1 as fol- 
lows, which is called the Dodgson's Determinant— Evaluation Rule (see 

det(A|i,„}) det(A) = det(An) det(A„„) - det(Ai„) det(A„i), (2.2) 

where A is an arbitrary matrix of order n and Aij is the minor of A by deleting the i-th 
row and the j-th column. 
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The following result shows the relation between the Pfaffian and the determinant. 



Lemma 2.3 (Godsil [5j) Let A be a square matrix of order n. Then 

A 
-A^ 



Pf 



det(v4). 



Let A = {ast)nxn be a skew symmetric matrix of order n and the corresponding 
directed graph. Suppose {i, j) is an arc in and hence aij > 0. Let be a directed graph 
with vertex set {1, 2, ... , n+1, n+2} obtained from by deleting the arc {i, j) and adding 
three arcs {i,n+ 1), (n + 1, n + 2) and (n + 2, j) with weights y/al[j, 1 and respectively 
(see Figures 1(a) and (b) for an illustration). For convenience, if aij = we also regard 

as a directed graph obtained from by adding three arcs {i, n + 1), (n+ 1, n + 2) and 
(n + 2, j) with weights 0, 1 and 0. The following lemma will play a key role in the proofs 
of our main results. 





(a) 



(b) 



Figure 1: (a) The directed graph G^. (b) The directed graph G^ 



Lemma 2.4 Suppose that A = (cisjnxn is a skew symmetric matrix and G*^ is the cor- 
responding directed graph. Let G^ be the directed graph with n + 2 vertices defined above 
and A{G^) the skew adjacency matrix of G^ . Then 

Pf{A) = Pf{A(G')). 

Proof Let G and G be the underlying graphs of G'^ and G^ , and let A(G'^) be the 
skew adjacency matrices of G"^. Hence A{G'^) = {ast)nxn and A{G^) = {bst){n+2)x{n+2), 
where 

ast ii 1 < s,t < n and {s,t) 7^ (i, j), (j,z), 

y/aij if (s,t) = (i,n + 1) or (n + 2, j), 
-^/a^ if (s, t) = (n + 1, i) or (j, n + 2), 
1 if (s,t) = (n + l,n + 2), 

-1 if (s,t) = (n + 2,n + l), 

otherwise. 
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By the definitions above, we have 



Pf{A)^Pf{A{G^)). 

Hence we only need to prove 

Note that, by the definition of the Pfaffian, we have 

Pf{A{G^))^ J2 Pf{A{G'))^ Yl ^' 

where M.{G) and M.{G) denote the sets of perfect matchings of G and G. 
We partition the sets of perfect matchings of G and G^ as follows: 

where Aii is the set of perfect matchings of G each of which contains edge e = {i.j), 
M.2 is the set of perfect matchings of G each of which does not contain edge e = 
A^i is the set of perfect matchings of G each of which contains both of edges + 1) 
and (n + 2,_7), and M.2 is the set of perfect matchings of G each of which contains edge 
(n + l,n + 2). 

Suppose TT is a perfect matching of G. If tt e then there exists uniquely a perfect 
matching n' of G — i — j such that tt = vr' U It is clear that there is a natural way 

to regard tt' as a matching of G. Define: vf = tt' U {{i, n + 1), {n + 2, j)}. Hence vf G M.i. 
Similarly, if tt e A42, we can define: 7f = tt U {{n + 1, n + 2)} and hence vf e M.2- It is 
not difficult to see that the mapping / : tt i — > W between M{G) and M{G) is bijective. 

Hence we only need to prove that for any perfect matching tt of G we have b^^ — hf^. By 
the definition of 7f, if tt = {(si,ti), (52,^2), • • • , j)) {si+i,U+i), . . . , (sf , t^)} G 

Ml, then 7f = {(si, ti), (s2, ^2), • • • , (^i-i, t^-i), (i, n+1), (n+2, j), (s^+i, ti+i), . . . , (s" , G 
Note that 

sgn{sxtx . . . si-iU-iij si+iti+i . . . sntn) = sgn{siti . . . si-iti-ii{n+l){n+2)isi+iti+i . • -Sfts), 

bsiti ■ ■ ■ bsi_iti_ibi{n+l)b(n+2)jbsi+iti+i ■ ■ ■ ^snbu — 

Qsitl ■ ■ ■ Qsj_iti_i-\/Ctij-\/0«jQs(+it(+i ■ ■ ■ '^Sntn — ^Sltl ■ ■ ■ Qs(_itj_i ^ij Qs;_|_it(_|_i ■ ■ ■ ^Sntn ■ 
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Thus we have showed that if tt G A^i then we have 6,^ = Similarly, we can prove that 
if TT e then we have 6^ = b^^. So we have proved that Pf{A(G'^)) = Pf{A(G^)), and 
the lemma follows. ■ 



3 New PfafRan identities 



We first need to introduce some notation. In this section, we assume that A = {aij)nxn 
is a skew symmetric matrix with n even. Suppose E = = 1, 2, . . . , /c} is a subset 

of [n] X [n] such that ii < ^2 < ■ ■ ■ < ^fe and i; < ji for 1 < / < /c. We define a new skew 
symmetric matrix E{A) of order n from A and E as follows: 



-^(^) = {bij)nxm bij 







if (i, j) ^ E and i < j, 
if (j, i) ^ and i > j, 
otherwise. 



By the definition of E{A), it is obtained from A by replacing all and {ji, i;)— entries 

with zeros and not changing the other entries and hence it is a skew symmetric matrix. 
For example, if ^4 = (ajj)4x4 is a skew symmetric matrix and E = {(1,4), (2,3), (3,4)}, 
then 

/ ai2 ai3 ^ 

-ai2 024 

-ai3 

-024 y 

Now, we can state one of our main results as follows. 



E{A) 



\ 



Theorem 3.1 Suppose A — {aij)nxn 'is a skew symmetric matrix of order n and E — 
= l,2,...,k} is a non empty subset of [n] x [n] such that ii < ^2 < ■ ■ ■ < 

ik, ii < ji for / G [k]. Then, for any fixed p G [k], we have 

Pf{E{A))Pf{A) = Pf{E,{A))Pf{E,{A))-r 
E ""nn [f(P, l)Pf{E{A\i^,j,})Pf{A{j^,i,}) - g(p, l)Pf{E(A){i^,i,y)Pf(A{j^j,y)] , 

l<l<k,l^p 

where Ep = E\{(ip,jp)}, Ep = {(ip,jp)},f(p,l) = s([n],ipji)s([n],jpii) and g(p,l) = 
s{[n],ipii)s{[n],jpji). 
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Proof Let be the corresponding directed graph of A defined as above, whose 
vertex set is [n]. Let be the directed graph with the vertex set [n + 2k] obtained 
from by replacing each arc between every pair of vertices ii and ji with three arcs 
{ii, n + 2/ — 1), (n + 2/ — 1, n + 21) and (n + 21, ji) with weights y/aij'i, 1 and ^/Oi^ if 
is an arc of G"^ and with three arcs {ji, n + 2l — l),{n + 2l — l,n + 21) and {n + 2/, ii) with 
weights y/cij^i, 1 and ^Ja^i if (j/, i;) is an arc of , respectively. For the case ai^j^ > 
for 1 < / < fc, Figure 2 (a) and (b) illustrate the procedure constructing G*^ from G^. 
Suppose A = A{G^) is the skew adjacency matrix of G^. 




(a) 



(b) 



Figure 2: (a) The directed graph G^. (b) The directed graph G^. 

Take a = [n], j3 = {n + l,n + 2, . . . ,n + 2k} = {xi\xi = n + i,l < i < 2k}. Take 
q = 2p — 1. Hence Xg = n + 2p — 1 and (— 1)'+'?+^ = (—1)'. By Lemma 2.2, we have 

2k 

Pf^{a)Pf^{aP) = Y,i~^)'PfA{oiX,xi)Ph{aP\x,xi). (3.1) 
1=1 



By the definitions of E{A) and G*^ and Lemma 2.4, we have 

PUa) = Pf{E{A)), Pf^{a(3)=Pf{A). 

We set 

au = -PfA{aXgX2l'-l)PfA{(yP\XgX2l'-l), 
hv = Pf-x{aXqX2l')PfA{(yP\XqX2l'), 

that is, 

a, = -Ph{a{n + 2p-l){n + 21' - l))P f^{aj3\{n + 2p-l){n + 21' - 1)), 

b, = Ph{a{n + 2p-l){n + 2l'))P f^{al3\{n + 2p-l){n + 21')). 
By Lemma 2.4, it is not difficult to see that 

h, = Pf^{aXgX,)PUa(3\{xgX,}) = P f{E,{A))P f(E;{A)). 
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(3.2) 



(3.3) 
(3.4) 

(3.5) 



Note that Op = 0. Hence we have 

2k 

1=1 

^Pf{E,{A))Pf(E;{A))+ Yl + (3.6) 

l<l'<k,l'^p 

Obviously, if ai^j^ = then theorem is trivial. Hence we may assume that Qi^j^ ^ 0. 

First, we prove that if Oj^jj, > then the theorem holds. From (3.2) and (3.6) it suffices 
to prove the following claim: 

Claim For any I' e [k] and I' ^ p, if ai^j^ > we have 

PfiEiAhp,,})PfiA{,,,,}). (3.7) 

Suppose ajpjp > 0. Then (ip, n + 2p — 1), (n + 2p — 1, n + 2p) and {n + 2p, jp) are three 
arcs of with weights -sjod^^i 1 a-nd ^fOd^- We need to consider two cases: 

(a) ai.ny > 0; 
(6) ai^,j^, < 0. 

If ciii/ji, > 0, then (i^/ , n + 2Z' — 1) , (n + 2Z' — 1, n + 21') and {n + 21' , ji/) are three arcs of 
G^ with weights y/cij^, 1 and y/cii^. Suppose X is a subset of the vertex set of G. Let 
be the directed subgraph of G^ induced by X and G[X] the underlying graph of 
GfX]^. Note that G[a{n + 2l' — l)(n + 2p— 1)] contains two pendant edges {ii',n + 2ii' — 1) 
and (ip, n + 2p — 1). Each perfect matching tt of G[a{n + 21' — l){n + 2p — 1)] can be 
denoted by tt = tt' U {{if, n + 2i;/ — 1), {ip,n+ 2p — 1)}, where tt' is a perfect matching of 
G[a\{ii/,ip}]. Set 

Pf{A(G[a{n + 2l'-l){n + 2p-l)Y))= _ E K: 

TTeMiG[a{n+2l' -l){n+2p-l)]) 

Pf{A{G[a\{z,,Zp}Y))= _E K', 

n'eM{G[a\{ii,,ip}]) 

where M.{G) is the set of perfect matchings of a graph G. By the definitions of 6^ and 
^tt', it is not difficult to see that 

= sgnip - l')s{[n],ipii').ya~a~b^', 

where sgn{x) denotes the sign of x. By the definition of E(A), we have 

Pf{A{G[a\{z,,Zp}r)) = Pf{E{A)^,^„^y). 
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Hence we have proved the following: 

Ph{a{n + 2l'-l){n + 2p-l)) = sgn{p - /')s([n], V«Wp«vJ,^/(^(^)fe,v})- (3-8) 
Similarly, we can prove the following: 
PUaP\{n + 2l' -l){n + 2p-l)) ^ sgn{p-l')s{[n],j,jv),/a-^ (3.9) 

Ph{a{n + 2l'){n + 2p - 1)) = sgn{l' - z^j,) V«vip«vi.'^/(^(^)feJ.}); (3-10) 

Ph{a(3\{n + 2l'){n + 2p- 1)) = sgn{l' - p)s(N, j>,OV%^p«vj-.,^/(A.p-v})- (3-11) 

Then (3.7) is immediate from (3.3), (3.4), (3.8) — (3.11). Hence ifa^,j^, > then the claim 
follows. 

If ai^,j^, < 0, then {ji',n + 21' — 1), (n + 2/' — 1, n + 21') and (n + 21', iii) are three arcs 
of with weights sJ—Od^n 1 and -^J—cii^,. Similarly, we can prove the following: 

PhHn + 21' - l){n + 2p- 1)) = syn(p - l')s{[nl Vr)^/^— a— P/(A(£;(A){,^,,-^,}); 

(3.12) 

P/^(Q;/3\(n + 2r-l)(n + 2p-l)) = s^n(p- r)s([n], j>,OV«^p.pS.v^/(Aw.'}); (3-13) 
P/^(a(n + 2/')(n + 2p - 1)) = s^n(/' - p)s([n], V«Wp%,'V^/(^(^)fe,v}); (3-14) 
Ph{aP\{n + 2r)(n + 2p - 1)) = s^n(r - p)s(N,ipjWA/W^^/(%pJ.})- (3-15) 

Then (3.7) is immediate from (3.3), (3.4), (3.12) — (3.15). Hence if a^^,j^, < then the 
claim follows. 

Hence we have proved that if ai^j^ > then the theorem holds. 

If ai^j^ < 0, we consider Pf{-A) and Pf{-E{A)). Note that {-A)i^j^ > 0. The 
result proved above implies that 

Pfi-EiA))Pfi-A) = Pf{E,{-A))Pf{E;{-A)) - a,^,^ J] (-a,,,Jx 

l<l<k,lj^p 

[f{p,l) X P/(E(-A){,^,,})P/((-A){,^,,j) -^(p,/)P/(E(-yl){,^,,j)P/((-yl){,^,,})] . 

(3.16) 

Note that by the definition of the Pfaffian we have Pf{-A) = (-l)tp/(^). By (3.16), 
we can show that we have 

Pf{E{A))Pf{A) = Pf{E,{A))Pf(E;{A))+ 
10 



l<l<k,l^p 

which imphes that if ai^j^ < then the theorem also holds. 

Hence we have proved the theorem. ■ 

Corollary 3.2 With the same notation as Theorem 3.1, for any fixed p G [k], 

Pf{E{A))Pf{A) = Pf{E,{A))Pf{E,{A))-r 

E [/(^^' - 9{V. l)Pf(E(A){j^,j,y)Pf(A{,^,,,})] . 

l<l<k,l^p 

Proof Let be the transpose of A. Note that Pf{A^') = The corol- 

lary follows immediately from Theorem 3.1 by considering the transpose of A. ■ 

The following result is a special case of Theorem 3.1 and Corollary 3.2. 

Corollary 3.3 Suppose A = {aij)nxn is a skew symmetric matrix of order n and E — 
{(^/) iOK = 1, 2, . . . , /c} is a non empty subset of [n] x [n] such that ii < ji < 12 < J2 < 
. . .ii < ji < . . . < ik < jk- Then 

Pf{E{A))Pf{A) - Pf{E,{A))Pf{E'M)) 

1=2 

= [Pf{E{A\,,,,})Pf{A{,,,,}) - Pf{E{A){,,,,})Pf{A{,,,^})\ . 

1=2 

Remark 3.2 The Pfaffian identities in Theorem 3.1 and Corollaries 3.2 and 3.3 express 
the product of Pfaffians of two skew symmetric matrices E[A) and A in terms of the 
Pfaffians of the minors of E{A) and A, where E{A) is a skew symmetric matrix obtained 
from A by replacing some non zero entries a^ji and Uj^ii of A with zeros. On the other 
hand, an obvious observation in the Pfaffian identities known before, which belong to the 
Pliicker relations, is that the related matrices are either a skew symmetric matric A or 
some minors of A. Hence the Pfaffian identities in Theorem 3.1 and Corollaries 3.2 and 
3.3 are completely new and different from the Pliicker relations. 

Example 3.1 Let A — (01^)4x4 and E = {(1,2), (3, 4)}. Then, by Corollary 3.3, we have 
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Pf 



Pf 



auaaiPf 




ai3 
-ai3 



4 Applications 



As applications of some results in Section 3, we obtain some determinant identities 
different from the Pliicker relations in Section 4.1 and we prove a quadratic relation for 
the number of perfect matchings of plane graphs in Section 4.2, which has a simpler form 
than the formula in |32J. 



4.1 New determinant identities 



by E[A] 



Jstjnxni 



where h 



St 



We first need to introduce some notation and terminology. Throughout this subsec- 
tion, we will assume A = {aij)nxn is an arbitrary matrix of order n and E = < 
^ ^ k} ^ [n] X [n], where ai^j^ ^ 0. Define a new matrix of order n from A and E, denoted 

ast if (s, t) ^ E, 

. In other words, E[A\ is an n x n 

otherwise. 

matrix obtained from A by replacing all entries ai^j^ for 1 < / < A; with zeros and not 
changing the other entries. For example, if A = (aij)4x4, E = {(1, 2), (2, 2), (3, 1)}, by the 
definition of E[A] we have 

/ an ai3 ai4 \ /an an a.iq a-iA \ 

0-21 a23 024 

a32 033 034 

y 041 0,42 ) 



E\A\ 



{(3, 4)} [A] 



ail Q'12 0-13 0-14 

O21 O22 O23 O24 

C^Sl ^^32 '^33 

y O41 O42 O43 O44 j 



Now, we start to prove the following: 
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, A . 

Lemma 4.5 If A — [aijjnxn is a matrix of order n and A* = I ^ | , then, for 

any i,j G [n], i ^ j, we have 

where Aij denotes the minor of A obtained by deleting the i—th row and j—th column 
from A. 

( B . ^ ^ 

Proof Note that A't- -i = , where -B is an (n — 2) x n matrix obtained 

( B . 

from A by deleting two rows indexed by i and J. Obviously, det I ^ I ~ 0- Hence 

2 [ b\ 

= det =0, which implies that 

[-Br o) 

Pf{Ayjy) — 0. Similarly, by Lemma 2.3 we can prove (ii) and (Hi). Hence the lemma 
follows. ■ 

Theorem 4.2 Let A = {aij)nxn be a matrix of order n and E = < I < k} a 

non empty subset of [n] x [n], where ii <i2 < ■ ■ ■ <ik- Then for a fixed p G [k] we have 

det{E[A]) det{A) 

^det{E,[A])det(E;[A]) - E det(£;[A],^,-J det(AdJ, 

l<l<k,l^p 

where Ep = E\{{ip,jp)} and Ep = {{ip,jp)}- 

A 

Proof Define: A* ^ \ ^ \ = (a*,.)2nx2n and E* = {{ii,n + ji)\l < I < k}. 



by Cayley Theorem we have PflAy^y) 



-A^ 



By Theorem 3.1, we have 

Pf{E*{A*))Pf{A*) = Pf{E;{A*))Pf{El{A*))+ 



l<l<k,lj^p '- 
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where f{p,l) = s{[2n],ip{n + ji))s{[2n], {n + jp)ii) and g{p,l) = s{[2n],ipii)s{[2n], {n + 
jp){n + ji))- It is not difficult to see that we have the following: 

<(n+,,) = <(n+n) = <^Hn. f(P, = -(-1)^^+^-^+^'+^'. (4.2) 

By Lemma 2.3 and the definitions of A* and E*[A*), we have 

Pf{E*{A*)) = (-l)^"Mdet(E[yl]), Pf{A*) = det(A), (4.3) 

Pf{E;{A*)) = (-l)^("-^)det(£;^[A]), Pf(E;{A*)) = (-l)^"("-^)det(:^[^]). (4.4) 
By (i) in Lemma 4.5, we have 

Pf(E*(A%^,,y) = 0, (4.5) 

and by (ii) and (iii) in Lemma 4.5, we have 

^/(^*(^*)fe,n+,,}) = (-l)§(-i)("-^) det(£;[A],^J, (4.6) 

^/(^{n+i.,.}) = (-l)^("-^)('^-^) det(A,,J. (4.7) 
The theorem is immediate from (4.1) — (4.7) and hence we have completed the proof of 
the theorem. ■ 

In the proof of Theorem 4.2, (4.1) is obtained from Theorem 3.1. Obviously, A corre- 
sponding identity to (4.1) can be obtained from Corollary 3.2. Similarly, by this identity 
we can prove the following: 

Theorem 4.3 Let A = {aij)nxn be a matrix of order n and E = < I < k} a 

non empty subset of [n] x [n], where ii <i2 < ■ ■ ■ ^ik- Then for a fixed p e [k] we have 

dei{E[A\)dei{A) 

= det(£;,[^])det(:^[yl]) - (-l)''+''+'^+'''«v..«iu.det(£;[A],,,gdet(^,,,J, 

l<«<fc,Z^p 

where Ep = E\{{ipjp)} and Ep = {{ipjp)}. 

The following result is immediate from Theorems 4.2 and 4.3. 

Corollciry 4.4 Let A — (ajj)„xn be a matrix of order n and E — {(i;,_7;)|l < I < k} a 
non empty subset of [n] x [n], where ii < ^2 < ■ ■ ■ < Then for a fixed p e [k] we have 

k 

^(-l)^'+^'a,,,-,{det(^[£;],,,J det(>l,,,J - det(£;[^],,,J det(A,,J} = 0. 

1=1 
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Example 4.2 Let A = (aij)3x3, E = {{1,1), (2, 2), (3, 3)} and p = 2. Then, by Theorems 
4-2 and 4-3, we have 



ail ai3 

021 0.22 0,23 
O3I O32 ^33 



ai2 ai3 
021 a23 
a-si 032 



ail 0-12 (213 

aai a23 

0-31 0,32 a33 



ai2 ai3 

^21 0-22 0-23 
agi 032 





021 


O23 




O12 


Ol3 




On 


Ol3 







012 


011(222 













~ 022033 














O3I 


O33 




O32 




021 


O23 




031 


032 





012 


Ol3 




O2I 


O23 




On 


012 







Ol3 


011O22 













— 022O33 














032 


O33 




O3I 




031 


032 




021 


O23 



4.2 Graphical edge— condensation for enumerating perfect match- 
ings 

Let M{G) denote the sum of weights of perfect matchings of a weighted graph G, 
where the weight of a perfect matching M of G is defined as the product of weights of edges 
in M. It is well known that computing M{G) of a graph G is an A^P-complete problem 
(see ^'J. Inspired by (2.2)-Dodgson's Determinant— Evaluation Rule, Propp 24 first 
found the method of graphical vertex-condensation for enumerating perfect matchings of 
plane bipartite graphs as follows: 

Proposition 4.2 (Propp \24\ ) Let G = {U, V) he a plane bipartite graph in which \U\ = 
\V\. Let vertices a, b, c and d form a 4— cycle face in G, a,c E U , and b,d E V . Then 

M{G)M{G - {a, b, c, d}) = M{G - {a, b})M{G - {c, d}) + M{G - {a, d})M{G - {b, c}). 
By a combinatorial method, Kuo ^1] generalized Propp's result above as follows. 



Proposition 4.3 (Kuo [14| ) Let G = {U,V) be a plane bipartite graph in which \U\ = 
\V\. Let vertices a, b, c, and d appear in a cyclic order on a face of G. 

(1) If a,c E U , and b,d E V , then 

M{G)M{G - {a, b, c, d}) = M{G - {a, b})M{G - {c, d}) + M{G - {a, d})M{G - {b, c}). 

(2) If a,b E U, and c,d E V , then 

M{G - {a, d})M{G - {b, c}) = M{G)M{G - {a, b, c, d}) + M{G - {a, c})M{G - {b, d}). 
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By Ciucu's Matching Factorization Theorem in Yan and Zhang [33j obtained 
a more general resuh than Kuo's for the method of graphical vertex-condensation for 
enumerating perfect matchings of plane bipartite graphs. Furthermore, Yan et al [32] 
proved the following results: 

Proposition 4.4 (Yan, Yeh and Zhang [32j ) Let G be a plane weighted graph with 
2n vertices. Let vertices ai, bi, 02, 62, • • • , Ofc, bk {2 < k < n) appear in a cyclic order on a 
face ofG, and let A = {ai, 02, ■ ■ ■ , a^}, B = {bi, b2, - ■ ■ , b^}- Then, for any j = 1,2, ■ ■ ■ , k, 
we have 

J2 M{G-aj-Y)M{G-A\{aj}-Y) = ^ M{G-W)M{G-A-W), 

YCB, \Y\ is odd wcB, \w\ is even 

where the first sum ranges over all odd subsets Y of B and the second sum ranges over 
all even subsets WofB,Y = B\Y and W = B\W . 

The following result, which is a special case of the above theorem, was first found by 
Kenyon and was sent to "Domino Forum" in an Email (for details, see |32j). 

Corollary 4.5 Let G be a plane graph with four vertices a, b, c and d (in the cyclic order) 
adjacent to a single face. Then 

M{G)M{G -a-b-c-d) + M{G -a- c)M{G -b-d) 

= M{G -a- b)M{G -c-d) + M{G -a- d)M{G -b-c). (4.8) 

By Ciucu's Matching Factorization Theorem in [2], Yan et al |S2] also obtained some 
results for the method of graphical edge-condensation for enumerating perfect matchings 
of plane graphs. In this subsection, by using the new Pfafiian identity in Corollary 3.3 we 
will prove a quadratic relation, which has a simpler form than the formula in [32], for the 
method of graphical edge-condensation for computing perfect matchings of plane graphs. 

We first need to introduce the Pfaffian method for enumerating perfect matchings 
IT^ . If G'^ is an orientation of a simple graph G and C is a cycle of even length, we 
say that C is oddly oriented in G^ if G contains odd number of edges that are directed in 
in the direction of each orientation of G. We say that is a Pfaffian orientation of 
G if every nice cycle of even length of G is oddly oriented in G'^ (a cycle C in G is nice if 
G — G has perfect matchings). It is well known that if a graph G contains no subdivision 
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of -^"3,3 then G has a Pfaffian orientation (see JE])- McCuaig McCuaig et al [20!, 
and Robertson et al. ||25| found a polynomial-time algorithm to show whether a bipartite 
graph has a Pfaffian orientation. 

Proposition 4.5 ( |12|, I17l |) Let he a Pfaffian orientation of a graph G. Then 

[M{G)f = det{A{G'), 
where A{G'^) is the skew adjacency matrix of G'^ . 

Remark 4.3 Let he a Pfaffian orientation of a graph G and A[G'^) the skew adjacency 
matrix of G^ . By Cayley Theorem and Proposition 4-5, we have 

M{G) = ±Pf{A{G^)), 

which implies that, for two arhitrary perfect matchings tti and 7C2 of G, hoth 67^1 (in'd h^^^ 
have the same sign. 

Proposition 4.6 (Kasteleyn's theorem, [111. I12|. I17| ) Every plane graph G has an 
orientation G^ such that every houndary face-except possihly the unhounded face has an 
odd numher of edges oriented clockwise. Furthermore, such an orientation is a Pfaffian 
orientation. 

Now we can prove the following result: 

Lemma 4.6 Let G he a plane graph with four vertices a,b,c and d (in the cyclic order) 
adjacent to the unhounded face. Let G^ he an arhitrary Pfaffian orientation satisfying the 
condition in Proposition 4-6 and A = A{G'^) the skew adjacency matrix of G'^. Then all 

Pf{A{aM4})PfiA)^PfiAa,b})Pf{A{c4}),PfiAa,c})Pf{A{b4}) andPf{A{a4})Pf{A{,^c}) 
have the same sign. 

Proof By (2.1) in Corollary 2.1, we have 

P/(A{a,6,c,4)^/(^)=^/(^{a,H)P/(A{,,,})-P/(A|„,,})P/(A{M})+^/(^{a,d})P/(A|5,4). 

(4.9) 

Obviously, ^{^,6,^,4, ^{a,6}, ^{c,d}, ^{a,c}, ^{6,4, ^{a,4 and are the skew adjacency 

matrices of G'^ — a — b — c — d, G'^ — a — b, G^ — c — d, G^ — a — c, G^ — b — d, G'^ — a — d 
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and C — b — c, respectively. Note that all the orientations G^ — a — b — c — d,G^ — a — 
b,G^ — c — d^C^ — a — c^C^ — b — d^C^ — a — d and G'^ — b — coiG — a — b — c — djG- 
tt — b.G — c — d,G — a — c,G — b — d.G ~ a — d and G — b — c satisfy the condition in 
Proposition 4.6 and hence are Pfaffian orientations. By Remark 4.3, we have 

M{G) = ±Pf{A), M{G-a-b-c-d) = ±Pf{A{,,,,,,a}). 
Hence we have proved the following: 

M{G)M{G -a-b-c-d)^ ±PfiA)Pf{A{aAc,d})- (4-10) 
Similarly, we can prove the following: 

M{G-a- b)M{G -c-d)^ ±Pf{A^a,b})Pf{A{c,d}y, (4-11) 

M{G-a- c)M{G - b - d) = ±P/(^{„,,})P/(A{,,,}); (4.12) 

M{G-a- d)M{G -b-c)^ ±Pf{A^a,d})PfiA{b,c}). (4.13) 

The lemma is immediate from (4.8) — (4.13). ■ 
Now we can start to state the main result in this subsection as follows. 

Theorem 4.4 Suppose G is a plane weighted graph with even number of vertices and 
the weight of every edge e in G is denoted by uje- Let ei = aifei,e2 = 0262, ■■■,6^ = 
Ofe&fe (k > 2) be k independent edges in the boundary of a face f of G, and let vertices 
Oi, fci, 02, b2, ■ ■ ■ , Ojfc, fcjfc appear in a cyclic order on f and let X — {ei\ i — 1, 2, ... , k}. 
Then, for any j — 1,2, . . . , k, 

M{G)M{G -X)^ M{G - ej)M{G - X\{ej})+ 

uje, ^ iJei[M{G-bj-ai)M{G-X -aj-bi)-M{G-bj-bi)M{G-X -aj-ai)]. 

Proof Note that Ci = Oi&i, 62 = 0262, ■ ■ ■ i^k — cik^k {k > 2) are k independent edges 
in the boundary of a face / of G. It suffices to prove the following: 

M{G)M{G -X)^ M{G - ei)M{G - X\{ei})+ 
ei J2 [M(G - 61 - ai)M(G -X-ai-bi)- M(G - 61 - bi)M(G - X - oi - a^)] , (4.14) 



k 
1=2 
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Since G is a plane graph, for an arbitrary face F of G there exists a planar embedding 
of G such that the face F is the unbounded one. Hence we may assume that vertices 
oi, 02, &2, ■ ■ ■ , cbk-i bk appear in a cyclic order on the unbounded face of G. Let T be a 
spanning trees containing k edges Cj's and let be an orientation of T such that the 
direction of each edge is from Oj to 6^ for i = 1, 2, . . . , k. Because each face of G can be 
obtained from T by adding an edge, it is not difficult to see that there exists an orientation 
G^ of G obtained from which satisfies the condition in Proposition 4.6. Hence all 
G'^ G^-X, G^-Cj, G^-X\{ej}, G'-ai^bj, G^-X-aj-bi, G^-b^-h and G^-X-ai-aj 
are Pfafiian orientations satisfying the condition in Proposition 4.6, the skew adjacency 
matrices of which are A, E{A), Ej{A), Ej{A), A{a^^b,}, E{A){aj,b,}, and E{A){ai,aj}, 

respectively, where E = {(ai,6j)|l < i < k},Ej = E\{ej} and Ej = E\Ej. Without 
loss of generality, we may assume that Oj = 2i — l,bi — 2i for i = 1,2, ... ,k, that is, 
E = {(1, 2), (3, 4), ... , (2A; - 1, 2A;)}. By Corollary 3.3, we have 

Pf{E{A))Pf{A) = Pf{EM))Pf{EM))+ 

k 

au J2 «2-i,2i [Pf{E{A){,,2i})Pf{A{2,2i-i}) - Pf{E{A){,,2i-i})Pf{A{2,2i})] . (4.15) 

i=2 

By a similar method to that in Lemma 4.6, we can prove that 

Pf{E{A))Pf{A) = ±M{G - X)M(G); (4.16) 

Pf{E^{A))Pf(E;{A)) = ±M(G - X\{e^})M{G - e,); (4.17) 

Pf{E{A)^,,2^})Pf{A{2,2^-l}) = ±M{G - X - - bi)M{G - b, - a,); (4.18) 

PfiEiA){^,2i-i})PfiA{2,2i}) = ±M{G -X-ai- ai)M{G -b,- h). (4.19) 

Since every perfect matching of G — X is also a perfect matching of G, by the definition 
of the Pfaffian, both Pf{A) and Pf{E{A)) have the same sign. Hence by (4.16) we have 

Pf{E{A))Pf{A) = M{G - X)M{G). (4.16') 

Similarly, we have 

Pf{E,{A))Pf(E[{A)) = M{G - X\{e,})M{G - e,). (4.17') 

Note that if vr' is a perfect matching of G — Oi — 6i — — fej [i ^ 1) then tt = vr' U 
{(ai, bi), {tti, bi)} is a perfect matching of G. By the definition of the Pfaffian, it is not 
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difficult to see that both 6,^ and b-,,' have the same sign, which imphes that both Pf{A) 
and Pf{A{a^^bi,a,M) have the same sign. Hence Pf{A)Pf{A{a^^b^,a,M}) ^ 0- Lemma 
4.6, we have 

P/(A|,,,,,l)P/(A{,,,,,}) > 0, P/(A|,,,,,^|)P/(A|,„,,}) > 0. (4.20) 

Since every perfect matching of G— X — ai — 6j is also a perfect matching of G—ai—bi, both 
Pf{E{A){aub,}) and Pf{A{ai,b,}) have the same sign. Similarly, both Pf{E{A){aua,}) and 
Pf{A{ai,ai}) have the same sign. Hence by (4.20) we have 

P/(E(A)|„,,,^})P/(%,,,,}) > 0, P/(E(A)|,,,,^})P/(%,.,,}) > 0. (4.21) 

From (4.18), (4.19) and (4.21), we have 

Pf{E{A){,,2^})Pf{A{2,2^-l}) = M{G - X - - hi)M{G -h,- a{)- (4.18') 

P/(E(A)|i,2.-i})P/(A|2,2.}) = M{G - X - ai - a,)M{G -h- h). (4.19') 

Note that ai2 = cUei and a2i-i,2i = ^e,- It is not difficult to see that (4.14) follows from 
(4.15) and (4.16') — (4.19'). Hence we have complete the proof of the theorem. ■ 

Remark 4.4 The formula in Theorem 4-4 for the method of graphical edge- condensation 
for enumerating perfect matchings of plane graphs has a simpler form than that in Theorem 
3.2 in J21 
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